Black hole formation in perfect fluid collapse 
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We construct here a special class of perfect fluid collapse models which generalizes the homoge- 
neous dust collapse solution in order to include non-zero pressures and inhomogeneities into evolu- 
tion. It is shown that a black hole is necessarily generated as end product of continued gravitational 
collapse, rather than a naked singularity. We examine the nature of the central singularity forming 
as a result of endless collapse and it is shown that no non-spacelike trajectories can escape from the 
central singularity. Our results provide some insights into how the dynamical collapse works, and 
into the possible formulations of the cosmic censorship hypothesis, which is as yet a major unsolved 
problem in black hole physics. 
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Black hole physics has attracted considerable attention in 
recent years and has witnessed rapid theoretical develop- 
ments as well as numerous astrophysical applications. It 
is to be noted, however, that while few exact static or sta- 
tionary models of black holes such as the Schwarzschild, 
Reissner-Nordstrom, and Kerr-Newman spacetimes are 
well-studied, the actual formation of black holes within 
the framework of a dynamical gravitational collapse pro- 
cess is not really an arena where much is known. 

In a realistic physical scenario, (stellar mass) black 
holes will be typically born when a massive star exhausts 
its nuclear fuel, and then collapses endlessly under the 
influence of its own gravitational field. Towards mod- 
eling such a physical process, a well-known model that 
has served as the basic paradigm in black hole physics 
is that of the Oppenheimer-Snyder spherically symmet- 
ric collapse solution where a dust cloud undergoes a 
continued collapse to form a black hole. Here the col- 
lapse initiates from a regular initial data, when there is 
no trapping of light (i.e. light rays from the star can es- 
cape to faraway observers). Subsequently, as the collapse 
advances the process of formation of an event horizon and 
closed trapped surfaces takes place, thus leading to the 
formation of black hole and the eventual spacetime singu- 
larity. The trapped surfaces and the event horizon form 
here well in advance to the epoch of the formation of the 
spacetime singularity, which is hence necessarily hidden 
within the black hole. 

Even though this collapsing Friedmann model already 
tells us a homogeneous dust collapse will always end in a 
black hole rather than a naked singularity, it should be 
noted that this scenario in fact has several limitations. 
For example, in this case the cloud has no pressures in- 
cluded, whereas any physically realistic collapse must in- 
clude pressures. Another restrictive assumption here is 
that the density profile is assumed to be strictly homo- 
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geneous in space, at all times through out the evolution 
of the cloud. For any isolated object such as a star, one 
may rather like to study a physically realistic density dis- 
tribution which would be typically higher at the center, 
and decreasing as we move away from the center. 

It is thus essential to study and analyze more general 
collapse situations in order to understand the black hole 
formation in more realistic stellar collapse scenarios. This 
is also essential in order to make any possible progress to- 
wards the cosmic censorship hypothesis 0], which broadly 
states that any physically realistic gravitational collapse 
must result into the development of a black hole. Such a 
conjecture has been absolutely fundamental to the theory 
of black holes, and has played a major role in astrophysi- 
cal applications of the black hole physics. This, however, 
remains a major unresolved open problem in general rel- 
ativity and black hole physics today. 

From such a perspective, we study here a specific class 
of collapse models where matter obeys the perfect fluid 
equation of state, and construct models where the col- 
lapse always necessarily ends in the formation of a black 
hole. The models we study here are somewhat special in 
that the mass function is assumed to be separable in the 
variables which are the physical radius of the cloud, and 
the time coordinate. However, this is a class which gen- 
eralizes the Oppenheimer-Snyder dust collapse models in 
two important respects, namely inhomogeneities of den- 
sity distribution are included and also non-zero pressures 
have been incorporated now. As the collapse always ends 
here in a black hole formation as we show, it is hoped that 
dynamical considerations such as these will provide some 
useful insights into physically realistic collapse and the 
actual process of black hole formation. It is not unlikely 
that it is only such dynamical considerations which would 
prove essential to resolve the issue of cosmic censorship. 
The model here may be of interest as it includes pressures 
which may be important in the later stages of collapse, 
and because the equation of state is that of a perfect 
fluid, which is physically a well-studied equation of state 
widely used in various astrophysical considerations. 

The fluid content of the cloud is in the form of a per- 
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feet fluid with an equation of state of the form p = kp, 
i.e. at all epochs the radial and tangential pressures are 
equal and isotropic, and are proportional to the density 
function of the cloud. Though the case of a general in- 
homogeneous dust collapse with k = can be completely 
solved 0, there are still a number of open questions re- 
garding the end state of a general perfect fluid collapse. 
Our purpose here is to examine a class of solutions of 
the Einstein equations for a spherically symmetric per- 
fect fluid to understand explicitly how an inhomogeneous 
density profile should behave in the later stages of col- 
lapse and near the singularity, so that the final state of 
the collapse would always be a black hole necessarily. 

The spacetime geometry within the spherically sym- 
metric collapsing cloud can be described by the metric in 
the comoving coordinates (t, r, 0, (f>) as given by, 



ds l 



z 2»(t,r) dt 2 + ^{t,r) dr 2 + R 2 ^ ^2 



(1) 



where dQ 2 is the line element on a two-sphere. The 
energy-momentum tensor for any matter fields of type 
I Q (this is a broad class which includes most of the 
physically reasonable matter fields, including dust, per- 
fect fluids, massless scalar fields and such others) is then 
given in a diagonal form, 



T 



-p(t, r); T; = p r (t, r); Tg = T* = pg(t, r) (2) 



The quantities p, p r and pg are the energy density, ra- 
dial and tangential pressures respectively of the cloud. 
We take the matter fields to satisfy the weak energy con- 
dition, i.e. the energy density measured by any local 
observer is non-negative. Then for any timelikc vector 
V\ we must have, 



T ik V % V* > 



which amounts to, 



P> 0; p+Pr > 0; p + p e > 



(3) 



(4) 



Now for the metric Q the Einstein equations take the 
form, in the units (8irG = c = 1) 



P = 



F' 
R 2 R' 



F 
R^R 



2(pe-p r ) R' 

P+Pr R 



-2R! + R^ + R^=0 
G H 



G- H = 1 

R 



(5) 
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where, 



The arbitrary function F — F(t, r) here has an in- 
terpretation of the mass function for the cloud, and it 
gives the total mass in a shell of comoving radius r on 
any spacelike slice t = const. We have F > from the 
energy conditions. In order to preserve the regularity 
at the initial epoch, we have F(U,0) — 0, that is the 
mass function should vanish at the center of the cloud. 
Since we are considering collapse, we have R < Q, i.e. 
the physical radius R of the cloud keeps decreasing and 
ultimately reaches R = 0. As seen from equation (6), 
there is a density singularity in the spacetime at R = 0, 
and at R' = 0. However, the later ones are due to shell- 
crossings and these weak singularities can be possibly re- 
moved from the spacetime J5j, so we shall consider here 
only the shell-focusing singularity at R = 0, which is 
the genuine physical singularity where all matter shells 
collapse to a zero physical radius. 

Now let us incorporate the perfect fluid form of matter, 
where the radial and tangential pressures are equal, and 
take the equation of state for the collapsing matter to be 



p r (t,r) =p e (t,r) = kp(t,r) 



(10) 



where k < 1 is a constant. Then the equations JSJ) and 
© become, 



F' 
R 2 R' 
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(ii) 



(12) 



Thus we see that there are five dynamical variables, 
namely p, iJj, v, R, and F, and there are five total 
field equations. Also, using the scaling independence we 
can write R(ti,r) — r at the initial epoch t = ti from 
where the collapse commences. The time t — t s (r) corre- 
sponds to the formation of the shell-focusing singularity 
at R = 0, where all the matter shells collapse to a van- 
ishing physical radius. 

Now let us assume that the mass function can be ex- 
plicitly written as a function of the physical radius R of 
the cloud and <, 



F(R,t) = R 3 M(R)Q(t) 



(13) 



G(t, r) = e- 2 ^(R') 2 ; H(t, r) = e- 2v {R) 2 (9) 



That is, we consider the class of mass functions F to be 
separable in R and t. Apart from that it is general in the 
sense that M is any C 2 function, whereas the function 
Q(t) yet to be determined by the field equations is a suit- 
ably differentiable function of t for t < t Sf> , where t = t S() 
is the time for the occurrence of the central singularity. 
One may consider equation (13) to be a somewhat strong 
assumption on the nature of the mass function. However, 
our basic purpose here is to construct a class of dynam- 
ical collapse models, where the two main constraints of 
the homogeneous dust collapse mentioned above are re- 
laxed, namely the density need not be homogeneous at 
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the initial epoch and later as well during the collapse 
evolution, and secondly we want to allow for non-zero 
pressures while considering a dynamical collapse situa- 
tion. As both these purposes are met by the above form 
of mass function as we have chosen here, it is adequate 
for our present purposes, as it allows us to construct ex- 
plicit collapse models which are more general, and which 
necessarily end up in a black hole as we shall see. 

Another requirement that is frequently imposed on 
physical grounds on the initial data from which the col- 
lapse evolves is that the physical variables such as the 
density and pressures are taken to be smooth or analytic 
functions at the initial surface. We can then write, 



M(R) 



1 1 

3 + 5 



M 2 R 2 



(14) 



Then from equation (|llf> we get p(r, t) = p(R, t) with, 

p(R, t) = (3M + RM iR )Q{t) = A(R)Q(t) (15) 

where the function A(R) is given by, 

A(R) = 1 + M 2 R 2 + ■■■ (16) 

As seen from above, at the initial epoch t = ti, the density 
function is given by, 



p(r,U) = p (r) = Q(U) [1 + M 2 r 2 + 



(17) 



Thus we see that the gradients of the density and pres- 
sures of the cloud vanish at the center at the initial epoch 
as required by the smoothness. Also, for the density to 
diverge at the singularity, we must have, 



lim Q(t) 



oo 



(18) 



It follows that for the given mass function the perfect 
fluid condition can be written as 



(k + l)QA+^MQ = 
R 



(19) 



The solution of the above equation determines the mass 
function completely. 

In order now to construct a class of collapsing solu- 
tions, let us consider the case when v = v(R), i.e. let the 
metric function v be a function of the physical radius R 
only. Again, this restriction is good enough for us as it 
allows us to construct the collapse models which include 
inhomogeneity and non-zero pressures, and which end 
up in black holes, generalizing the collapsing Friedmann 
models. A further useful feature of these choices may be 
considered to be that it allows to be incorporated a per- 
fect fluid with a reasonable equation of state, rather than 
any arbitrary forms of pressures (e.g. a purely tangential 
pressure, while assuming that the radial pressures iden- 
tically vanish) as is done some times when considering 
gravitational collapse. 

One can now integrate equation I) 1211 to get 



v(R) 



k 



k + 



-hx[CiA(R)} 



(20) 



Here C\ is a constant of integration. Now putting the 
value of H(t,r) in equation J7J and simplifying we get, 



R'G - 2Rv'G = 



(21) 



It is now possible to solve the above equation and the 
function G has the form. 



G(R) = A(R)~ 



(22) 



In other words, the above forms of p, v and G solve the 
Einstein's equations 0. Obtaining now the func- 
tion Q(t) will complete the solution. Putting in these 
functions in equation (|SJ) we get, 



R 



C 2 A—\J A— - l + R 2 MQ 



(23) 



where C 2 = C* +1 is another constant. The negative sign 
denotes the collapse condition R < 0. Finally, substitut- 
ing the values of the functions A and M we get, ignoring 
the higher order terms, as we are interested to find a 
solution close to the singularity, 



R 



C 2 R 



2kR 2 
k + 1 



1 



M 2 R? 



Q 



2kMo 



We note that in the equation 
changes sign at the value R 



(24), 



+ 1 
(24) 

the velocity R 
This would 



V(l + *0/2fc. 

correspond to a class of dynamic perfect fluid models 
where there is a bounce at the above value of the physical 
radius. However, since we are interested in the collapse 
models only presently, we do not consider this bouncing 
branch of the solutions. This can be achieved by fixing 
the boundary conditions suitably. For example, for the 
extreme value k = 1, this change corresponds to R = 1. 
Now as we have the scaling R — r at the initial epoch, 
this means that the boundary of the object r = r\, at 
the initial epoch is to be given by < < 1. Then at 
all later epochs this condition will be of course respected 
because R < 0, and the physical radius R monotonically 
decreases with t and will be less than one for all shells at 
all future times. For smaller values of fc, that is k < 1, we 
of course have larger values of the boundary of the cloud 
available, and in the extreme case k — 0, i.e. the dust 
collapse models, we can have arbitrarily large rj, without 
the velocity ever changing sign, or the cloud can be as big 
as we want, and there will be no bounce at all possible 
in the dust case. 

Now let us solve the equation 119|) close to the space- 
time singularity at R = 0. In this approximation, 
A(R) — ► 1 and M(R) — ► |. Using these approximations 
and equation (|24l) in equation 119fl we get 



(k+l)Q- 



1 



3CW¥ 



= 



2fc 
fe+1 



M 2 



(25) 



Considering that M 2 < and solving the above equation 
with the boundary condition Q(t So ) — > oo as pointed out 
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earlier, we get 
Q(t) = -a+ 
where, 



2a 



[exp{-^C 2 (fc + l)a(t SQ - t)} - 1] 

(26) 



6k\M 2 



1 



(27) 



Thus we see that the above Q(t) is solution to Ein- 
stein's equations in the vicinity of the singularity with 
respect to the given forms of p, v and G. Now we can 
also solve for the metric function R, which is the physical 
radius for the cloud. Using equation (|2~lf> we get, 



R(r,t) = f(r)e- B ^ 



(28) 



where f(r) is an arbitrary function of r. To avoid any 
shell crossing singularity we consider / to be an increas- 
ing function of r and since the area radius of the geomet- 
rical center of the cloud vanishes we must have /(0) = 0. 
The function B(t) is given as, 



B(t) = C 2 d% 



2k 
k + 1 



M 2 dt 



(29) 



As noted earlier, the spacetime singularity occurs at 
R = 0. We now need to decide if the singularity in 
the present case is necessarily covered within an event 
horizon of gravity (which is the case of a black hole for- 
mation), or it could be visible to faraway observers in 
the spacetime. The way to decide this is to examine if 
there are any future directed families of null geodesies 
which go out to external observers in future, and which 
in the past terminate at the singularity. If such families 
do exist, then the singularity is naked, which in princi- 
ple can communicate with outside world, and in the case 
otherwise we have a black hole forming as the end state 
of collapse. We thus need to consider the existence or 
otherwise of such families of paths from the singularity. 

With the form of R as given above in equation (28), 
and as R < 0, the singularity happens at a time t = t s 



when the physical radius for all the shells with different 
values of the comoving coordinate r becomes zero. In 
other words, there is a simultaneous collapse of all the 
shells to singularity, and as t — > t So all the shells labeled 
by the coordinate r collapse simultaneously to the singu- 
larity at R = 0. This necessarily gives rise to a covered 
central singularity at R = 0, r = 0, as there are no out- 
going future directed non-spacelike geodesies coming out 
from the same. Because, if there were any such outgoing 
geodesies, given by say t = t(r) in the (t, r) plane, which 
came out from t = t s ,r = 0, then the time coordinate 
must increase along these paths, which is impossible as 
there is complete collapse at t = t s , and there is no space- 
time beyond that. Hence no values t > t s are allowed 
within the spacetime which does not extend beyond the 
singularity. Thus, the collapse gives rise necessarily to a 
black hole in the spacetime. 



Our main purpose here has been to generalize the 
homogeneous dust collapse scenario to include non-zero 
pressures, and the inhomogeneities of density and pres- 
sures, which is physically more realistic situation, while 
ensuring that the collapse end state is a black hole only. 
While the perfect fluid collapse models we considered 
here allow for inhomogeneities in density and pressure 
profiles, and do necessarily give rise to black holes as we 
have shown here, it should be kept in mind, as we have 
noted above, that the classes of mass functions and the 
velocity profiles for the collapsing shells as determined 
by the choice of metric function v{R) considered here are 
rather special. It is an open problem to explore if we 
could generalize these assumptions further, and if so to 
what extent, and still continue to get black holes only 
and not the naked singularities as the final end product 
of gravitational collapse. The point is it is known, for ex- 
ample, for inhomogeneous dust collapse 0] that as long as 
the inhomogeneities are within certain limits, the result 
of collapse is a black hole. However, beyond that critical- 
ity of inhomogeneities, the collapse could end in a naked 
singularity. Investigating further specific, but physically 
more realistic models, may illustrate better such features 
of gravitationally collapsing configurations. 
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